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Abstract: The purposes and possibilities of future diffractive structure function 
measurements at HERA are presented. A review of the current range and accu- 
racy of the measurement of F^ 3 \(3, x p , Q 2 ) is presented and an estimate of the 
precision of future measurements is given. A feasibility study is performed on the 



Tj - 1 measurement of the structure functions F^ )( - 4: \p,x lp ,Q 2 ,t), F® charm , R D< ^\ R D ^ 

and Ff. Included in this study are estimates of the integrated luminosity required, 
the analysis techniques to be employed and values of systematic error that could be 
expected. 

a: 

1 Introduction 

The first pioneering measurements of the deep inelastic diffractive cross section at HERA have 
^ ■ yielded great insight into the mechanism of diffractive exchange and the partonic structure of 
the pomeron ^, H, The precision of present measurements is, however, inadequate for the 
study of many quantitative aspects of this picture. The origin of the gluonic structure of the 
pomeron and the interface between soft and hard physics are unresolved questions that would 
benefit from the higher luminosity offered by HERA in the years to come. In addition, the 
substantially different partonic structure of the pomeron in comparison to usual hadrons makes 
diffraction a challenging test for perturbative QCD ||. Furthermore, it has been suggested 
that the emerging data from HERA may provide fundamental insight into the non-perturbative 
aspects of QCD, leading to a complete understanding of hadronic interactions in terms of a 
Reggeon Field Theory 0. 

In this paper we study the measurements that can be made of inclusive diffractive cross 
sections with the luminosities which may be achieved in the future running of HERA. We 
evaluate what precision is possible for these measurements, and attempt to highlight those 
measurements which will be of particular relevance to achieving significant progress in our 
understanding of QCD . 



x To appear in the proceedings of the workshop "Future Physics at HERA", DESY, Hamburg, 1996. 



2 Experimental Signatures of Diffractive Dynamics 



A prerequisite for any precision measurement of a hadronic cross section is that it must be de- 
fined purely in terms of physically observable properties of the outgoing particles. Since there 
is, as yet, only the beginnings of a microscopic understanding of the mechanism responsible for 
diffractive dynamics, and since there are a plethora of wildly different predictions for how the 
diffractive cross section should depend on all of the kinematic variables, it is the experimental- 
ist's job to provide a well defined measurement unfettered by any ad hoc assumptions about 
how the cross section should behave. 

For the interaction ep —>■ eX (X = Anything) the cross section is usually measured differen- 
tially as a function of the two kinematic variable x and Q 2 defined as follows: 
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where q and P are the 4-momenta of the exchanged virtual boson and incident proton respec- 
tively, and y is the fractional energy loss of the electron in the rest frame of the proton. The 



differential cross section 
F L = F 2 - R/(l + R) by: 



dxdQ 2 



may then be related to the proton structure functions F 2 and 
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Such an interaction may be further characterised by dividing the total hadronic final state 
(i.e. the final state system excluding the scattered lepton) into two systems, X and Y, separated 
by the largest gap in true rapidity distribution of particles in the 7*-p centre of mass system, as 
shown in figure |]. When the masses of these two systems, Mx and My, are both much smaller 




Figure 1: Schematic illustration of the method of selecting events to define a diffractive cross-section, showing 
the definitions of the kinematic variables discussed in the text. The systems X and Y are separated by the 
largest gap in true rapidity in the 7*-p centre of mass system. The system Y is nearest to the direction of the 
incident proton in this frame of reference. 

than the 7*p centre of mass energy W then a large rapidity gap between the two systems is 
kinematically inevitable, and the interactions are likely to be dominantly diffractive. 

The two additional kinematic degrees of freedom may be specified by introducing the vari- 
ables t, (3 and x p : 

t=(P-Y) 2 , (3) 
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The hadronic cross section to be measured is therefore Ax 1qT^ y Av The presence of a rapidity 
gap may then be identified experimentally by use of the calorimeters and forward detectors of 
the experiments. For sufficiently small values of Mx, the system X will be completely contained 
in the central calorimeters, allowing an accurate determination of (3 and x p . Tagging of energy 
in the very forward direction allows both the mass of system Y and the magnitude of \t\ to be 
constrained, but does not allow a direct measurement of t. 

A leading proton spectrometer (LPS) "tags" particles which follow closely the path of the 
outgoing proton beam separated by only a very small angle. By tagging the particle at several 
points along it's trajectory and with detailed knowledge of the magnetic field optics between 
the interaction point and the tagging devices, it is possible to reconstruct the original vector 
momentum of the particle. Thus it is possible to select a sample of events containing a "leading" 
proton with E p ~ E^ am . An LPS permits an accurate measurement of t, and provides an 
independent method of selecting interactions sensitive to diffractive dynamics. Since there is 
no need for all the hadrons of system X to be contained within the main calorimeter it is 
possible to make measurements at higher values of x p , so that subleading contributions to the 
cross section may be confirmed and further investigated. 

Regge phenomenology has been highly successful in correlating many features of the data 
for both elastic and inelastic hadron-hadron cross sections J7|. At high energies (low x p ) the 
data are dominated by the contribution from the leading trajectory of intercept a(0) > 1. The 
most simple prediction of Regge phenomenology is that at fixed Mx the cross section should 
vary oc x~ p n where n is related to the intercept of the leading trajectory by n = 2a(t) — 1.0 The 
size of the rapidity gap is kinematically related to the value of x p by Ar] ~ — lna^ ||. It is 
easy to show that at fixed M x the rapidity gap distribution is then related to a{t) by: 

(LV oc (e- A ") 2 - 2 "«. (6) 
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Hence for a(t) < 1 the production of large rapidity gaps is exponentially suppressed, whilst for 
a{t) > 1 it is exponentially enhanced 0. With x p fixed, then extensions of the simple Regge 
model are necessary to predict the Mx dependence. For sufficiently large masses, it is perhaps 
possible to use Muller's extension of the optical theorem |IJ ("Triple Regge") to achieve this 



aim, and for small masses the vector meson dominance model |TT[ may be appropriate. 

There is much evidence from hadron-hadron data [0, and some from ep data || that ad- 
ditional ("sub-leading") trajectories are necessary to reproduce the measured hadronic cross 
sections. The Regge formalism permits t dependent interference between some of the differ- 
ent contributions rendering problematic any theoretical predictions for the dependence of the 
measured hadronic cross sections on M x , xjp, W 2 or t. 

Estimates of the diffractive contribution to both deep-inelastic and photoproduction cross 
sections have been made using the M x dependence at fixed W M. 



2 It is worth pointing out that at Q 2 > then x F is the basic Regge variable, not W 2 as is sometimes 
assumed. 



3 Future Measurements of F 2 



The 1994 HI data 0, |J, corresponding to ~ 2 pb -1 , have allowed a measurement of F 2 to 
be made in the kinematic region 2.5 < Q 2 < 65 GeV 2 , 0.01 < (3 < 0.9 and 0.0001 < x F < 0.05. 
For Q 2 < 8.5 GeV 2 this was achieved by taking data with the nominal interaction point shifted 
by +70 cm allowing lower angles to be covered by the backward electromagnetic calorimeter. 

The dependence of on ip was found not to depend on Q 2 but to depend on /3, 

demonstrating that the factorisation of the 7*p cross section into a universal diffractive flux 
(depending only on xp) and a structure function (depending only on (3 and Q 2 ) is not tenable. 
These deviations from factorisation were demonstrated to be consistent with an interpretation 
in which two individually factorisable components contribute to F^ 3 \ These two components 
could be identified with pomeron (IP) and meson contributions oc x p np ,x~ p nM where rip = 
1.29 ± 0.03 ± 0.07 and n M = 0.3 ± 0.3 ± 0.6. Scaling violations, positive with increasing logQ 2 
for all (3, were observed and could be interpreted in terms of a large gluon component in the 
diffractive exchange, concentrated near x 9 /jp = 1 at Q 2 ~ 2.5 GeV 2 [[12], |13|, |[ |]. 

Given the significant progress in the understanding of diffractive dynamics that has been 
achieved with the existing data, the goal of future measurements is twofold: to extend the kine- 
matic regime of existing measurement, and to achieve highest possible precision, particularly 
where existing measurements have uncovered interesting physics. In the 1994 measurement, 
with 5 bins per decade in both Q 2 and x p , and 7 bins in f3 between 0.01 and 0.9, then in the 
interval 8.5 < Q 2 < 65 GeV 2 there were an average of ~ 100 events per bin, corresponding to a 
statistical accuracy of ~ 10%[]. The different sources of systematic error for the 1994 measure- 
ment are shown in table [I], along with an estimate of the level to which they may be reduced in 
the future. The largest single error arises from the combination of limited Monte Carlo statis- 
tics (10%) and different possible final state topologies leading to varying corrections for finite 
efficiency and resolution (10%). The latter contribution was estimated from the difference in 
the correction factors calculated for two possible 7*iP interaction mechanisms: 

• a quark parton model process in which the 7* couples directly to a massless quark in the 
exchange with zero transverse momentum, 

• a boson-gluon fusion process in which the 7* couples to a gluon in the diffractive exchange 
through a quark box. 

Numerous experimental measurements of diffractive final state topologies at HERA are now 
available || which constrain the data to lie between these two possibilities. Therefore, it is 
reasonable to assume that the error arising from a lack of knowledge of the final state topologies 
may be reduced by a factor of ~ 2 such that it no longer dominates the total error. Monte 
Carlo statistics can obviously be increased. Therefore, it is reasonable that in the future a 
measurement may be made with a total systematic uncertainty of 5% or less. To reduce the 
statistical error to 5% (3%) would require 8pb _1 (22 pb" 1 ) of data. 

Less luminosity is required to achieve the same statistical precision at lower Q 2 . However, 
to reach the lowest possible Q 2 with the widest possible range in x (and hence in x p ) it 

3 The variation of the cross section with Q 2 of approximately Q~ 4 was partially offset by an increasing bin 
size with increasing Q 2 . 



HI 1994 Preliminary 


HI Future 


Error Source 






Main Calo' Hadronic E Scale: ±5% 


3% 


< 1% 


Backward Calo' Hadronic E Scale: ±15% 


3% 


< 0.5% 


Backward Calo' Elec. E Scale: ±1.5% 


5% 


<3% 


Tracking Momentum Scale ±3% 


2% 


< 1% 


Scattered Lepton Angle: ±1 mrad 


2% 


1% 


t dependence: e 6±2i 


1% 


0.5% 


xp dependence: x"^' 2 


3% 


1% 


(3 dependence 


3% 


1% 


Mx/xp resolution 


4% 


2% 


Background (photoproduction and non-ep) 


0.1% 


< 0.1% 


MC Statistics/Model Dependence 


14% 


4% 



Table 1: Sources of systematic error for the HI 1994 Preliminary measurement of F® . 
The results of calculations to estimate the extent to which these uncertainties may be reduced 
are shown in the right hand column. These calculations take cognisance of the new SPACAL 



backward calorimeter installed by HI [14], and rely upon future measurements made using the 



LPS and forward neutron calorimeter based upon a minimum luminosity of 10 pb . 

is advantageous to run with the interaction vertex shifted forwards along the proton beam 
direction. We calculate that 2 pb -1 of such data would give a statistical accuracy of 5% in 
the region of overlap between nominal and shifted vertex data, allowing a precise cross check 
of the two analyses. It is important to note that the HERA magnets are not optimised for a 
shifted vertex configuration, and so such data take twice as long to collect as for the nominal 
configuration. 

Theoretically, and consequently experimentally, the region of high (3 is of particular interest. 
The evolution of with Q 2 at high (3 is expected to depend crucially upon which evolution 
equations are pertinent to diffractive dynamics. In particular, a DGLAP |15|] QCD analysis JT^, 



H f|] demonstrates that the HI data are consistent with a large gluon distribution, concentrated 
near x g /jp = 1 at Q 2 ~ 2.5 GeV 2 . In this case then F®^ should begin to fall with increasing 
Q 2 at (3 = 0.9 for Q 2 3> 10 GeV 2 . The presence of a sufficiently large "direct" term in the 
evolution equations would lead to a indefinite increase with Q 2 [fL6|1 . Thus a measurement of 
F®^ at high (3 to the highest possible Q 2 is desirable. 

At fixed ep centre of mass energy y/s the range of (3 that may be accessed decreases with 
increasing Q 2 such that: 

(3 > ^ . (7) 

' arr rnax„ max v ' 
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The acceptance for a measurement of F®^ at f3 = 0.9, based upon the interval 0.8 < (3 < 1, 
therefore extends to a maximum Q 2 of ~ 3500 GeV 2 for xp < 0.05, where the cross section is 
likely to be dominantly diffractive. To achieve 10% statistical precision in a measurement in 
the interval 2000 < Q 2 < 3500 GeV 2 (bin centre Q 2 = 2600 GeV 2 ) would require ~ 200 pb" 1 . 
However, a measurement of 30% statistical precision, requiring only ~ 20 pb -1 , would already 
be significant in light of theoretical calculations which differ by more than a factor of 2 in this 
region. 



4 Measurement of F 2 



For particles with E <C Ep (here Ep is the proton beam energy) then the HERA magnets 
separate them from the proton beam allowing particles to be tagged for t ~ t m i n . For particles 
with E ~ Ep then only those particles with transverse momentum Pp > 0.07 GeV 2 may be 
tagged ||18|| . Consequently, for diffractive measurements then the acceptance in t is limited to 
t ^ —Pp/(1 — x F ) = —0.07 GeV 2 . For a process with a highly peripheral t dependence then this 
limitation results in the loss of 30 to 40% of the total cross section. The geometrical acceptance 
of the current ZEUS detectors is ~ 6% |0, giving an overall acceptance in the region x p < 0.05 



of ~ 4%. 

Tantalising first measurements of the t dependence of deep-inelastic events with a leading 



proton have been presented by the ZEUS collaboration ||18|| . The observed dependence on t 
of da/dt oc e bt , b = 5.9 ± 1.2^^7 GeV -2 , lends strong support to a diffractive interpretation of 
such interactions. The measurements of b differential in x p , (3 and Q 2 that will be possible 
with increased luminosity are eagerly awaited. 

A preliminary measurement F®^ using an LPS has also been made In order to use 



measurements of F®^ to estimate = ft max dt F^^ it is necessary to extrapolate from the 

measured region at high \t\ into the region of no acceptance to t m i n . To make this extrapolation 
with the minimal assumptions about the t dependence, it is necessary to have at least three bins 
in t for each bin in /3, Q 2 and x p . This means that ~ 150 pb -1 of data are required to make a 
measurement of similar accuracy to the existing data for F®^ ■ Obviously it would be possible 
to study the t dependence in aggregate volumes of phase space and make a measurements 
with a factor 3 fewer statistics, relying upon model assumptions about the variation of the t 
dependence with j3, Q 2 and x p . Even with 150 pb -1 of data, the extrapolation into the t region 
of no acceptance relies on the assumption that the t dependence is the same in the measured 
and unmeasured regions. It is not clear that the resultant theoretical uncertainty will be less 
than the 5% uncertainty in the overall normalisation resulting from proton dissociation for 
existing measurements of F^ 3 \ 

The primary purpose of the LPS is the diffractive "holy grail" of measuring F®^ (/?, Q 2 ,x p ,t). 
Measurements of any possible "shrinkage" of the forward diffractive scattering amplitude (in- 
creasing peripherality of the t dependence with decreasing x p ) are likely to have unique power 
for discriminating between different theoretical models of diffractive dynamics |5|. In addition, 
the ability to select subsamples of events in which there is an additional hard scale at the 
proton vertex is of great theoretical interest ||. It is therefore of the utmost importance that 
the 100 to 150 pb -1 of data necessary to make an accurate measurement of F^ i \/3,Q 2 ,x p ,t) 
are collected by both collaborations with the LPS devices fully installed. 

The LPS also allows measurements at higher x p for which the rapidity gap between the 
photon and proton remnant systems X and Y becomes too small to observe in the main 
calorimeter. This will allow the contributions to the measured hadronic cross section from 
subleading trajectories to be further investigated. Combining information from the scattered 
lepton and proton will allow the invariant mass Mx to be reconstructed in the region x p ^> 0.05, 
where it would otherwise be impossible. 

Tagged particles with E -C E h p am will provide information about the way in which protons 
dissociate into higher mass systems. Particles thus produced are kinematically unlikely to gen- 



erate a significant background in the region x p < 0.05, but for larger x p then this background 
will become increasingly important. 

5 On the Determination of FF 



A fundamental question in the study of hard diffractive processes is the extent to which 
perturbative QCD dynamics may be factorised from the proton vertex. Some recent calculations 
of diffractive cross sections f^, ^TJ consider the interaction in terms of a hard (perturbative) 
phase producing a coloured partonic system which subsequently interacts with the colour field 
of the proton via the exchange of non-perturbative gluons. Such "soft colour interactions" [21 
allow colour to be exchanged between the photon and proton remnant systems such that in 
some fraction of the events both will become colour singlet states separated by a potentially 
large rapidity gap. In such models the evolution of the effective parton distribution functions 
(PDFs) will be driven by the evolution of the PDFs of the proton. Alternatively, it is possible 
that the presence of a large rapidity gap will confine any evolution dynamics to within the 
photon remnant system X. The effective PDFs will then depend only upon f3 and Q 2 , and not 
upon x p as in the former case. 

The 1994 HI data have been demonstrated to be consistent with a factorisable approach 
in which a large gluon distribution is attributed to the pomeron in the region x p < 0.01 or 
over the whole kinematic range for the sum of two individually factorisable components 
At next to leading order (NLO) then a large gluon distribution G JP (/5, Q 2 ) necessitates a large 
longitudinal structure function, F F (P,Q 2 ): 
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A prediction for F L ([3, Q 2 ) based upon a NLO QCD analysis of the data could be tested directly 
since the wide range in x p accessible at HERA leads to a wide variation in the eP centre of 
mass energy ^/s e jp = ^x p s. This means that in factorisable models, at fixed (3 and fixed Q 2 the 
same partonic structure of the pomeron may be probed at different values of x p , corresponding 
to different values of y: 
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If the dependence of the diffractive structure function can be factorised such that 
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and F F is non-zero, then a measurement of 
value of F^ 3 ^ by the factor 8(y, /3, Q 2 ): 



assuming F^ =0 is lower than the correct 



F 2 D (3) {Measured) 2(1 - y) + y 2 /[l + R F (J3, Q 2 )} 
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Thus under the assumption of factorisation, a measurement of F F (j3,Q 2 ) is possible by mea- 
suring the extent of the deviation from a simple x~ n extrapolation from higher x p . The size 
of this effect is shown for different values of R F = F F / (F F — F F ) in figure 0. Also shown 
is the range covered by the existing 1994 data and the range that could be covered by future 
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Figure 2: The expected fractional deviation S(x p ) of from an extrapolation from high 

x F for R F = (horizontal solid line), R F = 0.5 (dashed line), R F = 1.0 (dotted line) and 
R F = 5.0 (dash-dotted line) for Q 2 = 10 GeV 2 , (3 = 0.5. The vertical solid line indicates the 
kinematic limit of y — 1. The larger box indicates the area covered by the preliminary 1994 
HI measurements, which extended up to y ~ 0.5 with an accuracy of ~ 20%. The smaller box 
represents a future measurement with a total error of 7% extending to y = 0.8. 

measurements at HERA. It is clear that in order to see this effect it is necessary to make a 
measurement at y values above 0.7 with a total error of less than 10%. For high values of f3 the 
extrapolation of into the high y region of interest should not be significantly compromised 
by the possible presence of subleading trajectories, which contribute only at high x p (low y). 

A comparison between the values of F F determined from such apparent deviations from 
factorisation and the values expected from a QCD analysis of F®^ constitutes a powerful test 
of both the validity of factorisation and the applicability of NLO QCD to diffraction at high 
Q 2 . 

6 Measurements of R D ^ and R D ^ 



Since an evaluation of F F using the techniques described in the previous section require 
theoretical assumptions concerning factorisation, such an analysis is clearly no substitute for 
a direct measurement of the ratio of the longitudinal to transverse diffract ive cross sections, 
R D (x p , j3). A good measurement of this quantity is vital for a full understanding of the diffrac- 
tive mechanism and should provide an exciting testing ground for QCD. There is at present no 
theoretical consensus on what values to expect for R D , although all models suggest a substan- 
tial dependence on (3 with most suggesting an extreme rise as /3 — *■ 1 [H]. A measurement of 
R D to any precision leads us into unexplored territory. 

Measurements of R D have so far been restricted to DIS exclusive vector mesons produc- 
tion |22] by a direct measurement of the polarisation of the final state resonance. This method 
could perhaps be used for the bulk data if the directions of the final state partons could be 
inferred, but is likely to be difficult due to the problems of running jet analyses on low mass 
final states. Instead we investigate a slightly modified version of the method used to determine 



R(x,Q 2 ) for inclusive DIS [23| 



The general form relating the structure functions F 2 D and to the ep differential diffractive 



cross section can be written in analogy to the inclusive cross sections [24 
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where R D ^ = (F 2 D{4) - 2a;F 1 D(4) )/(2a:F 1 D{4) ) = a%/a%. Although a measurement of R D ^ as 
a function of all 4 variables is the most desirable measurement and must be an experimental 
goal, statistical limitations are likely to mean that initial measurements must be made without 
a reliance on a leading proton spectrometer (LPS) and so t is not measured. In this case we 
define F 2 D(3) and R D ^ as 



d3< - 4TO2 ;i-y + n/1 , on £ „, J F^ {x ,Qi, Xr) , (13) 



dx r dxdQ 2 xQ 4K y 2(l + R D ^(x,Q 2 ,x 1P )) 

In this case R D ^ is the ratio of the longitudinal to transverse cross section only if R D ^ has 
no dependence on t. 

Analysis of equation 13 reveals that in order to make a measurement of R D independent of 
at least two ep cross sections must be compared at the same values of x, Q 2 and x F but 
different values of y. This is achieved by varying the ep centre of mass energy, y/s. There are 
of course many possible running scenarios for which either or both beam energies are changed 



to a variety of possible values. A full discussion on this point is given in [23]. For the present 



study we examine the case when the proton beam is roughly halved in energy from 820 GeV to 
500 GeV and the electron beam remains at a constant energy of 27.5 GeV so that data is taken 
at the 2 values of s of s = 90200 GeV 2 and s = 55000 GeV 2 . This setup allows for a reasonable 
luminosity at the low proton beam energy and enables systematic uncertainties concerned with 
detection of the scattered electron to be minimised. In this scheme we make a measurement 
of the ratio of the ep differential cross sections, r = cr^/a^, for two values of y, yu and yi 
(corresponding to the high and low values of s) for fixed x, Q 2 , x p and (if measuring R D ^) t. 
Equation [T2| or [13] is then used to determine R D . 

It is also apparent from equation O that in order to have the greatest sensitivity to R D 



measurements must be made at the highest yi Q possible (and thus lowest electron energy). This 
is illustrated in figure |, where it can be seen that for values of yi Q = 0.5 (or lower) there is 
little change of r for different values of R D . The upper limit in yi Q is crucially dependent on 
the ability of the detectors to discriminate and veto against photoproduction events in which 
a pion is misidentified as an electron. Experience has shown, however, that for the diffractive 
events the low mass of the final state reduces the chance of faking electron when compared to 
the more energetic non- diffractive events. For this study we take a central value of yi Q = 0.8 
with a lower (upper) bin limit of 0.75 (0.85) so that good electron identification for energies 
above 4.15 GeV is assumed. 

The kinematic range of the measurement projected onto the x-Q 2 plane is shown in figure [| 
for both CMS energies. To ensure that the scattered electrons are well contained within the 
backward detectors we place a maximum 6 e cut of 174°. This restricts us to Q 2 > 5 GeV 2 
and x > 10~ 4 . In order to ensure good acceptance in the forward region we impose a cut of 
x p < 0.01. 

For low electron energies the kinematic variables are well reconstructed and have good 



resolutions if the electron only method is used [25, 26 . Since the major problem with this 




Figure 3: Dependence of the ratio, r, of the ep cross sections at s — 90200 and s = 55000 with 
R D for various values of y at s = 55000. 




Figure 4: A projection of the kinematic range onto the x-Q 2 plane for a) a proton beam of 
E p = 500 GeV and b) a proton beam of E p = 820 GeV for an electron beam of 27.5 GeV. The 
shaded region represents the proposed region of study. Also shown is the restriction on the 
kinematic range imposed by a maximum 9 e cut of 174°. 



measurement will be in reducing the statistical error we select bins as large as possible whilst 
maintaining enough bins to investigate any variation of R D with the kinematic quantities. A 
suitable choice would be 4 bins per decade in x, 4 bins in (3 and if the LPS is used 2 bins in t. 
The bins in (3 and t are optimised so as to contain approximately equal numbers of events at 
each x value. Identical bins in these variables are used for both CMS energies. 
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logio x 


Number of Events 
without LPS 


Number of Events 
with LPS 


E P = 820 GeV 


E P = 500 GeV 


E P = 820 GeV 


E P = 500 GeV 


-4.125 


41 


28 


0.82 


0.56 


-3.875 


36 


25 


0.72 


0.50 


-3.625 


19 


13 


0.38 


0.27 


-3.375 


9 


6 


0.18 


0.12 



Table 2: The estimated number of events in each bin for an integrated luminosity of 1 pb~ 
for the 2 proton beam energies and an electron beam energy of 27.5 GeV, assuming 4 bins per 
decade in x, 4 bins in (3 and (for measurements with the LPS) 2 bins in t. R D was set to 0.5. 



In order to estimate the statistical errors on the measurement we used the RAPGAP gen- 
erator [271 with a fit to the measured HI F®^ 1 P| at s — 90200 GeV 2 and used equation IT2 



to determine the statistics at s — 55000 GeV 2 . We assumed 100% efficiency for measurements 
made without any LPS and 4% efficiency for those with. The expected number of events per 
integrated luminosity in each bin is summarised in table |2] for an example R D = 0.5. 

For systematic errors we estimate an error of 8{r)/r of 5%. This error is conservatively 
evaluated by taking the estimated error on F^ (see above) and assuming any improvement that 
arises from taking a ratio is offset by increased uncertainty in the photoproduction background 
and radiative corrections. 

An example of the sort of precision which may be obtained for a measurement of R D ^ and 
^>r>(4) j g sriown j n figme |5|. For this study we assumed that many more data would be obtained 
at the high CMS energy. It can be seen that for an integrated luminosity of 10 pb _1 at the 
lower s value a measurement of R D ^ is statistically dominated with an error around 60% if 
R D ( 3 ) = 0.5 for the lowest value of x. For an integrated luminosity of 50 pb _1 at the lower s 
value statistical and systematic errors become comparable and R D ^ can be measured to 40% 
accuracy. For measurements of R D ^ very high integrated luminosities are required - at least 
a factor of 50 is needed for a similar precision to R D ^. 

7 Measuring Ff charm 



Since the leading mechanism in QCD for the production of charm quarks is the boson gluon 
fusion process, the diffractive charm structure function F^ charm is very sensitive to the gluonic 
component of the diffractive exchange. It is important to establish whether the measured 
F® c harm * s consistent with that expected from a QCD analysis of the scaling violations in 
F^ 3 \ In addition, it has already been observed in the photoproduction of J/tp mesons that 
the charm quark mass provides a sufficiently large scale to generate the onset of hard QCD 
dynamics. The extend to which the charm component of F®^ 1 exhibits a different energy 
{x p ) dependence to that of the total will provide insight into the fundamental dynamics of 
diffraction. 

The method used here for tagging charm events uses the D* + decayQ D* + — > D°7i^ low — ► 
(K~7r + )7rf low . The tight kinematic constraint imposed by the small difference between the D* + 

4 Charge conjugate states are henceforth implicitly included. 
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Figure 5: The estimated errors for an example central value of R D = 0.5 for a) 10(500) pb at 
s = 55000 GeV 2 and 50(2500) pb" 1 at s = 90200 GeV 2 and b) 50(2500) pb" 1 at s = 55000 GeV 2 
and 250(12500) pb" 1 at s = 90200 GeV 2 for a measurement of R D{3) (R D ^). The inner error 
bar represents the statistical error and the outer the statistical and systematic error added in 
quadrature. 

and D° masses means that the mass difference AM = M(KinT s i 0W ) —M(Ktt) is better resolved 
than the individual masses, and the narrow peak in AM provides a clear signature for D* + 
production. The chosen D° decay mode is the easiest to use because it involves only charged 
tracks and because the low multiplicity means that the combinatorial background is small and 
that the inefficiency of the tracker does not cause a major problem. 

A prediction of the observed number of events is obtained using RAPGAP with a hard 
gluon dominated pomeron structure function taken from a QCD analysis of 0. The 

cross section predicted by this model for D* ± production in diffractive DIS is compatible with 
the value measured in j4|. The acceptance of the detector is simulated by applying cuts on 
the generated direction (9) and transverse momentum (p±) of the decay products and on the 
energy of the scattered lepton (E' e ). The p± cut used is 150 MeV, which is approximately 
the value at which the HI central and forward trackers reach full efficiency. This cut has 
a major influence on the acceptance for D* + mesons, because the momentum of the slow 
pion tt s i ow is strongly correlated with that of the D* + , so a D* + with p± much less than 
150 MeVx M£)*+/M 7r + « 2 GeV cannot be detected. The p_i_-dependence of the acceptance is 
shown in figure |6|a. There is no obvious way of extending the tracker acceptance to lower p±, 
so this cut is not varied. 

Figure |6|b shows the average acceptance for a D* + over the region 10 GeV 2 < Q 2 < 100 GeV 2 
and all values of f3 and for p± > 2 GeV. It can be seen that extending the angular coverage 
from the present 25 < 9 < 155° range in the backward direction to 170° in conjunction with 
lowering the scattered lepton energy cut used in present analyses significantly improves the 
acceptance, especially at low x p . Figure [5] shows the number of D*^ which one might expect 
to detect in the low- and high-Q 2 regions with a total integrated luminosity of 750 pb" 1 . It can 
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Figure 6: The acceptance for D* + in the region 10 GeV 2 < Q 2 < 100 GeV 2 , shown function 
of (a) p± and (b) x p . The continuous line in (b) shows the results with central tracking only 
and a requirement E' e > 10 GeV. The other lines show the effect of extending track coverage 
in the backward direction and including E' e down to 5 GeV. 



be seen that even with this large integrated luminosity, cross section measurements can only 
be made with an accuracy of 10% in this binning. 
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Figure 7: The number of D expected to be observed in 750 pb in the range (a) 10 GeV" < 
Q 2 < 25 GeV 2 and (b) 50 GeV 2 < Q 2 < 100 GeV 2 , predicted using RAPGAP with a hard- 
gluon-dominated pomeron. 



Whilst one waits for 750 pb 1 , it may be worthwhile attempting to increase the statistics 
by investigating other decay modes which are experimentally more demanding. The D* + decay 



to D°ir^ low has a branching fraction of nearly 70% |53| and is the only decay mode giving a 
charged track in addition to the D° decay products. However, the D° decay to K~tt + has a 
branching fraction of slightly less than 4% ||28|| , so there is clearly room for a large improvement 
in statistics if other channels can be used. For example the use of a silicon vertex detector 
close to the interaction point, such as that already partly installed in HI, should enable the 
secondary vertex from the decay of the D°, which has a decay length cr = 124 fim [p8| , to be 
tagged. This could be used to extend the analysis to other channels, including semileptonic 
decays, which are otherwise difficult to use. The gain in statistics, neglecting the vertex-tagging 
inefficiency, can be up to a factor of ~ 20 (if all channels are used), with a further factor of ~ 2 
available if inclusive D° production is used rather than relying on the D* + decay. 



8 Summary and Conclusions 



We have argued that a precise hadron level definition of the cross section to be measured is 
essential in order that future high precision measurements of diffractive structure functions may 
be treated in a consistent theoretical manner. Although 20 pb" 1 of integrated luminosity will be 
enough to achieve precision in the measurement of at moderate Q 2 , in excess of 100pb _1 
is necessary to make a comprehensive survey of diffractive structure within the kinematic limits 
of HERA. An attempt to determine will be an important tool in establishing the validity of 
both factorisation and NLO QCD in diffractive interactions. A direct measurement of R D ^ is 
shown to be feasible with 10 pb _1 of integrated luminosity taken at lower proton beam energy. 
A substantial integrated luminosity is demonstrated to be necessary to complete an exhaustive 
study of the diffractive production of open charm, although statistics can be markedly improved 
by exploiting the full range of possible charm decays. 
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